Two-spin relaxation of P-dimers in Silicon 
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We study two-electron singlet-triplet relaxation of donor-bound electrons in Silicon. Hyperfine 
interaction of the electrons with the phosphorus (P) nuclei, in combination with the electron-phonon 
interaction, lead to relaxation of the triplet states. Within the Heitler-London and effective mass 
approximations, we calculate the triplet relaxation rates in the presence of an applied magnetic field. 
This relaxation mechanism affects the resonance peaks in current Electron Spin Resonance (ESR) 
experiments on P-dimers. Moreover, the estimated time scales for the spin decay put an upper 
bound on the gate pulses needed to perform fault-tolerant two-qubit operations in donor-spin-based 
quantum computers (QCs). 

PACS numbers: 03.67.Lx,72.25.Rb,76.20.+q, 76.30.Da 



Solid state based quantum information processing is a 
rapidly developing field, bridging quantum information 
science and condensed matter physics. Various schemes 
have been proposed to use spins as quantum bits (qubit), 
with the most prominent examples using electron spins in 
quantum dots 1 and 31 P nuclear spins in Si. 2 Significant 
and exciting experimental progress have since been made 
to demonstrate coherent manipulation and measurement 
of spins in semiconductor nanostructures£~— 

It has been long demonstrated experimentally that sin- 
gle donor electron and nuclear spins in bulk Si have ex- 
tremely long coherence and relaxation times ^ Recent ex- 
periments on isolated donors show an electron spin de- 
coherence time T2 of 0.6 s at 1.8 K and 0.3 Tesla mag- 
netic field (up to 10 s at press tvme)Z& This coherence 
time scale could be modified for donors near an interface, 
although measurements for single donors near an oxide 
interface have obtained a promising donor electron spin 
relaxation time of T\ ~ 6 s at a magnetic field of 1.5 
Tesla£ 

The study of spin coherence properties of phosphorus 
dimers (two shallow P donors that are close to each other 
so that they are exchange coupled) in Si is relatively lim- 
ited compared to those for single donors, because bulk 
experiments are ensemble averages over many different 
dimer configurations and couplings, while two-donor arti- 
ficial devices are still in their infancy. P-dimers are an im- 
portant ingredient in the spin-QC architectures that use 
the nuclear or electron spins of 31 P in Si as qubits^^ In 
these proposals exchange interaction between two bound 
electrons lead to an effective coupling between the cor- 
responding electron or nuclear spins, which is used to 
implement two-qubit operations. Therefore, understand- 
ing spin coherence of P-dimers is crucial for the study of 
spin based quantum computing in Si:P. 

Here we calculate singlet-triplet relaxation rates of 
donor-bound electrons in bulk silicon within the Heitler- 
London approximation^ In the absence of pure dephas- 
ing (such as that due to 29 Si nuclear spins), such relax- 
ation forms one of the decoherence channels for two-spin 
states of a dimer, and would give an upper bound to two- 
spin coherence. This relaxation process could also lead 



to broadening of resonance peaks in an ESR measure- 
ment of P-dimers, making it an observable effect within 
the currently available spin resonance technology. 

Theoretical Model. We consider two electrons which 
are bound to two phosphorus donors (P dimer), sepa- 
rated by a distance 2d, as shown in Fig. [TJi. The crys- 
tallographic axes of silicon are denoted by x, y and z. 
The direction of the applied magnetic field is along Z 
(not necessarily z) axis, which defines the quantization 
axis in the spin space. The Hamiltonian of the coupled 
electrons- nuclei system is then given by H = Hq + Vhf 
where 
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g (g n ) is the effective g- factor of the electron (nucleus), 
Mb (aO is the electron (nuclear) magnetic moment, J de- 
notes the Coulomb exchange interaction between the two 
electrons, and Ahf =0.2 //eV is the hyperfine coupling 
constant between a phosphorus nucleus and its bound 
electron? 2 * The superscripts e (n) refers to the electron 
(P nucleus), and 1 and 2 label the first and the second 
electron (nucleus). Notice that the total electron spin is 
not a good quantum number due to the hyperfine interac- 
tion Vhf- For example, if the two nuclear spins are anti- 
aligned, the two electrons would experience an inhomo- 
geneous magnetic field, which mixes the electron singlet 
and triplet states^ This mixing then allows a phonon- 
mediated relaxation between the two-electron spin states, 
and is at the core of the present study. 

Using the Heitler-London approximation and averag- 
ing over the electron orbitals, we simplify the above 
Hamiltonian by introducing the averaged hyperfine cou- 
pling A: (lP as \VHF\i>as) ^ (tp s \V HF \ips) * \A{a ln + 
<T 2 ™) • (<X le + <X 2e ) and {lp aS \V H F\lp s ) = (lps\VHF\lpas) - 

\A[cr ln — <T 2n )-(er le — <r 2e ), where the symmetric and an- 
tisymmetric electronic orbital wave functions are denoted 
by ip s and tp a s > respectively. Considering only the ground 
orbital state on each donor, ip a ,as — {^(j"!)^^) ± 
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V'L( T *2)'i/'-R( T 'i)}/'V / 2, where i/jl (iI>r) is the ground orbital 
wave function of the left (right) donor electron to be 
specified later. We note that the first excited state of 
the donor bound electron is about 10 meV above the 
ground states*, and we do not consider dimers that are 
too tightly coupled (d > 8 nm), so that the Heitler- 
London approximation is justified for our calculation. 

Within the above mentioned approximations, the 16 x 
16 spin Hamiltonian of the two electrons and two nu- 
clei reduces to six clusters (one quartet, five doublets, 
and two singlets) that are block-diagonal. The effec- 
tive Hamiltonian for the quartet cluster of the electrons 
Hilbert space, in {S e S n , TfiTtf, T+T™, T^T™} bases, 
reads 
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FIG. 1: (a) Schematic of a P-dimer and the corresponding 
(hydrogenic) electrons envelope functions, (b) Energy levels 
of the two electrons bound to a P-dimer in the quartet clus- 
ter, see Eq. @- Mixing of the donor electron eigenstates 
is achieved via the hyperfine coupling to the P nuclear spins, 
which causes the triplet-singlet relaxation of the electron spins 
via emission of acoustic phonons. Here V refers to different 
relaxation rates from triplet states and 4J is the splitting be- 
tween T n e and S e . 



Ji = J + 2gn B B + 2g n fi n B - 2A, (4) 
J 2 = J - 2gfi B B - 2g n p n B - 2A, (5) 

where {S,T } = -± \U)T lit) and {T + ,T_} = {|ft 
), 144-)} are the singlet and the triplet spin states of the 
electrons and the phosphorus nuclei, as shown in Fig. \I]p. 
The effective Hamiltonians for the remaining doublets are 
given by 

# (2) = 2 /),{5 e T ",T ^™}, 
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The fully polarized states, {T£T£} and {T^T™}, are de- 
coupled from the remaining states. 

In the following we neglect the Zeeman splitting of the 
phosphorus nuclear spins because it is about three orders 
of magnitude smaller than the electron Zeeman splitting. 
We also assume that the averaged hyperfine coupling is 
much smaller than the exchange energy A -C J to keep 
our perturbation treatment valid, so that level alignment 
is qualitatively given by Fig. [T]b. 

The electron-phonon interaction. Direct relaxation of 
electron spin states are generally forbidden since electron- 
phonon interaction conserves spin. However, as shown in 
Eq. (J3J, electron spin singlet and triplets are mixed via 
the hyperfine interaction in a P-dimer. At low energies 
(J < 1 meV), which is our main concern here, only acous- 
tic phonons contribute to the spin relaxation processes. 
Moreover, only deformation potential coupling, consist- 
ing of a dilation part and a shear part, contribute to the 



relaxation processes in Si (or Ge). The Hamiltonian for 
an electron in one of the valley minima of Si i o 13 ' 14 

Hg v =E d Tr{e} + E u (k-e-k), (6) 

where e is the strain tensor due to phonons and k is the 
unit vector along the direction of one of the [100] con- 
duction band minima in the reciprocal space. E u is the 
shear deformation potential due to uniaxial strain along 
the [100] direction and Ed + E u is the volume deformation 
potential 14 . The electron-phonon interaction for a donor 
electron in the ground state is then given by-i£ 

H% = (E d + \e u ) J2 q J—^ K + « t -, J ) ; (?) 

where a q annihilates a phonon with momentum q, p = 
2330 kg/m 3 is the silicon density, ui q — viq is the en- 
ergy of the longitudinal acoustic phonon, vi is its corre- 
sponding velocity and we only consider the longitudinal 
phonons which have the dominant contribution to the 
electron-phonon interaction. We note that in Si, Ed ~ 5.0 
eV and E u = 8.77 eV-i 4 - 

The relaxation rate. To calculate two-spin relaxation, 
we need the two-electron orbital wave functions for the 
P-dimer. The conduction band of Si has six degenerate 
minima located close to (but not at) the edge of the first 
Brillouin zone ^ In the effective mass approximation, the 
ground state wave function of the donor bound electron 
is given by the homogeneous superposition of these six 
valleys 

*( r ) = iE^w^w. (8) 

Mr) = u„(r)e ik ^, (9) 

where /i = {±x, ±y, ±z} is the valley index, <f>u{r) are the 
Bloch wave functions and F^(r) are their corresponding 
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TABLE I: Relaxation time Ti (ms) for T e T n -» S e S n (see 
Fig. [T]b) as a function of J (meV) and d (A) for substitutional 
donors located along [001], [Oil] and [111] axes, ao = 5.43 A 
is the lattice constant in Si, ai = ao/v2 and 02 = v3ao/4. 
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envelope functions. To find a simple analytical expression 
for the spin decay rate, we first assume Gaussian envelope 
functions for the donor electrons 
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where a and b are the effective Bohr radii, and the re- 
maining envelope functions are obtained by cyclic change 
of x, y and z in Eq. ([TO]) . Using the Fermi golden rule and 
Eq. <[7|) , we calculate the relaxation rate (the one-phonon 
process) from Tjf to S e and obtain 

T To ^s « a 1 {3 d + ^S u fA 2 J(n qo +l)\I\ l , (11) 

1 _ . 
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where a\ ~ 100 a 2 — 25/ir 2 ph 4 vf , q = AJ/hvi, I — 
(\f , L(r)|5'fl(r)) is the overlap integral of the right and 
left donor wave functions, n q is the Bose-Einstein dis- 
tribution of phonons, and T G is the relaxation rate for a 
Gaussian envelope function. The relaxation rates Tt + ->s 
and Tt_^s have the same form as in Eq. (1111) but with 
the following replacement J —> (Jx,2 + 3J)/4 [see Eqs. 
(I4I5[) ]. Clearly, V has an explicit linear dependence on 
the exchange energy J and a quadratic dependence on 
the hyperfine coupling A. 

Since in general J oc I 2 , one might conclude that T 
should scale as J 3 . However, the exchange integrals 
are non-trivial functions of the distance between the 
donor sites d and have oscillatory dependence on donor 
positions^ so that a simple polynomial scaling is gener- 
ally not available. 

A more realistic envelope function F^r) for a donor 
should be hydrogenic with anisotropy. For example, 
Kohn and Luttinger proposed a variational form 12 
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FIG. 2: (Color online) ^most-cubic dependence of the relax- 
ation rate T of the symmetric two-electron triplet state Tq on 
the singlet-triplet splitting 4J, for substitutional donors along 
[001] direction, see Eq. (|11|) . An external magnetic field B — 1 
T is applied, in order to split the triplet states. Note that for 
TI state, we took only those Js which were larger than the 
Zeeman energy. 



With this envelope function there is no closed analytical 
form for the overlap integrals and the electron-phonon 
matrix elements. We therefore calculate the relaxation 
rates numerically. Table I shows our results for the re- 
laxation rate of the unpolarized triplet to singlet state, 
where we have taken into account the oscillatory behav- 
ior of J as a function of the donors separation 2d. The 
relaxation times change from tens of /is to hundreds of 
ms as we vary the distance between the substitutional 
donors 2d (and their relative orientation) from 8 to 12 
nm. We also note that the value of the exchange, as a 
function of the donor sites, strongly depends on along 
which symmetry axes of the crystal the two donors are 
aligned. le We find an almost cubic dependence of the re- 
laxation rate on 4J for donors along [001] direction, as 
shown in Fig. [21 On the other hand, there is no reli- 
able polynomial fit for the other two directions [011] and 
[111]. 

Our calculation for J is based on the plane wave 
approximation^- of the Bloch wave function <p(r) — 
Uk(r)e lk r , where we only keep the plane wave part and 
neglect the lattice periodic part Uk(r). The numerical 
values of J do not differ significantly (specifically in the 
Heitler-London approximation) from the more elaborate 
approaches which incorporate the full Bloch structure of 
the Kohn-Luttinger electron wave function o 17 i 18 . There- 
fore, our final results for the relaxation times T% give the 
right order of magnitude due to the linear dependence 
of the rates on J. We stress that our calculations are 
valid in the regime of relatively large exchange energies 
A <C J ■ In other words, although J is oscillatory as a 
function of d, it remains always finite (over a large range 
of values) and larger than the hyperfine coupling in [100] 
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FIG. 3: (Color online) Two-spin To —¥ S relaxation times as 
a lunction of the exchange splitting for donor pairs, which are 
located about 10 nm apart from each other and along different 
symmetry axis. The position of one of the donors varies within 
its nearest neighboring sites, in a radius 5 Angstroms. The 
solid and dotted lines show that the resulting data is within 
two polynomial curves of J 3 and J. 



direction. In contrast, for donors located apart along 
[110] or [111] axes, there are certain values of d which 
result in (almost) zero exchange energy^ where our per- 
turbation scheme breaks down. 

In reality, it is challenging to implant donors at spe- 
cific lattice sites. Figure |3] shows the two-spin relax- 
ation time, calculated with Eq. (jlip . as a function of 
exchange splitting J, for donor pairs which are about 10 
nm apart (along different symmetry axis), with the po- 
sition of one of the donors varying within a radius of 5 
Angstrom. Our results show that even with this small 
fluctuation in the donor location, there could be a three- 
orders-of-magnitude variation in both exchange splitting 
and triplet-singlet relaxation rate, ranging from 0.1 ms 
to 100 ms. 

The main result of our calculations is that even in puri- 
fied samples (with no 29 silicon isotopes), the two-electron 
states can relax via the combined effect of the hypcrfinc 
interaction between phosphorus nuclei and the electrons, 
and the electron-phonon interaction^ Furthermore, this 
relaxation mechanism is strongly dependent on the donor 
separation and their orientation along different symme- 
try axes of the crystal. Our results provide an upper 
limit to the duration of the two-qubit gates in the Kane 
architecture^ where the two phosphorus nuclei are en- 
tangled indirectly via the bound electrons through the 
hyperfine and exchange interaction. Due to the strong 
dependence of the relaxation time of the electronic triplet 
states on the location (and orientation) of the donors, 



the physical realization of donor-based schemes requires 
a careful control over the donor positions. 

The data presented in Table I and Figs. (|2l3p . as 
well as the expressions for the relaxation rate F, allow 
the identification of optimal separations and orientations 
for substitutional donors which mitigate this relaxation 
mechanism, yet still permit sufficiently fast two-qubit 
operations. The key is that the relaxation rate is pro- 
portional to J\I\ , while the speed of two-qubit gates 
depends on J linearly, whether for electron or nuclear 
spin qubits. The additional dependence on inter-donor 
overlap means that at larger inter-donor separations, or 
for those pair-positions that have reduced overlap due to 
valley interference, a larger ratio of J/T should be avail- 
able, so that fault-tolerant two-qubit gates are possible. 
For nuclear spin qubits more specifically, we note that 
the use of the electron spins as intermediaries for two- 
qubit gates 2 could have significant impact on gate fidelity. 
While isolated P nuclear spins are known to have out- 
standing coherence properties^ the electron spin relax- 
ation mechanism studied here could put a much stronger 
constraint on the gate speed, making the electron spin co- 
herence properties the ultimate gauge to determine the 
feasibility of nuclear spin qubits for quantum information 
processing. 

In the electron spin spectroscopy for P dimers, transi- 
tions between the electron spin triplet states are detected, 
while singlet states are not directly involved. The finite 
relaxation calculated here from triplet states to the sin- 
glet state presents a leakage for the triplet populations, 
thus should lead to a broadening of the ESR signal for 
the dimers. 

We note that the effects of the spin-orbit interaction 
can be safely neglected due to the small spin-orbit cou- 
pling in silicon^ However, for shallow donors located 
close to an interface, the interface roughness could po- 
tentially lead to a sufficiently strong electric field near 
the bound electrons, which would in turn induce a non- 
negligible extrinsic spin-orbit interaction. Due to the lack 
of a comprehensive model of the interface and its effect 
on the spin of the electrons, we did not cover this issue 
in this work and instead focused on bulk properties. 

In conclusion, we have calculated the relaxation rates 
of the electronic triplet states to the singlet state in 
P-dimers. These rates have nontrivial dependence on 
the singlet-triplet exchange splitting, and are strongly 
anisotropic. This relaxation mechanism can be studied 
in P-dimers spin spectroscopy experiments by measur- 
ing the linewidth of the ESR peaks. Finally, our results 
can be easily extended to similar materials (like Ge) with 
different valley degeneracies. 

We thank financial support from NSA/LPS through 
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